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ABSTRACT 



1 INTRODUCTION 



Collisional de-excitation rates of partially deuterated molecules are different from the 
fully hydrogenated species because of lowering of symmetry. We compute the collisional 
(de)excitation rates of ND2H by ground state para-H2, extending the previous results for He- 
lium. We describe the changes in the potential energy surface of NH3- H2 involved by the pres- 
ence of two deuterium nuclei. Cross sections are calculated within the full close-coupling ap- 
proach and augmented with coupled-state calculations. Collisional rate coefficients are given 
between 5 and 35 K, a range of temperatures which is relevant to cold interstellar conditions. 
We find that the collisional rates of ND2H by H2 are about one order of magnitude higher 
than those obtained with Helium as perturber. These results are essential to radiative transfer 
modelling and will allow to interpret the millimeter and submillimeter detections of ND2H 
with better constraints than previously. 

2 POTENTIAL ENERGY SURFACE 



Doubly deuterated ammonia has been detected for the first time by 
[RouefF et ar] ( |2000^ towards the dark cloud L 134N via its 1 10 - loi 
ortho and para transitions at 1 10 GHz and subsequently in the pro- 
tostellar environment L1689N by jLoinard et al.H2001[ l. The sub- 
millimeter fundamental transitions have then been detected at the 
Caltech Sub-millimeter Observatory by |Rouetf et al. 



pOSjl; 



Lis 



Gerin et al. 



|et al.| ( |200 6^ and on the APEX antenna in Chile by 
( |2006^ towards the Barnard 1 molecular cloud and L1689N. ND2H 
was shown to be a sensitive tracer of the physical conditions of this 
star forming region. In a previous paper l ,Machin & Roueff|2007^ 
some of us have presented the collisional formalism and com- 
puted the collisional excitation cross-sections and rate coefficients 
of ND2H by He by introducing the appropriate changes in the po- 
tential energy surface of NH3 - He, computed by Hodges & Wheat- 
|ley| ( |200l| . This new accurate intermolecular potential energy sur- 
face (PES) for NH3-He has also been used previously to reevalu- 
ate the collisional excitation of NH3 and NH2D by He ( Machin &| 



[Roue fi^ 2005 , Scifoni et al. 2007 , M achin & Roueff|2006t . Here, we 
present results for ND2H in collision with para-H2. H2 is the main 
constituent of dark clouds environment where ND2H has been de- 
tected. Collisional excitation is dominated by interaction with low 
temperature H2; this motivated our present study, aimed particu- 
larly at the low temperature excitation schemes of ND2H. 

The paper is organized as follows. In section|2]we briefly de- 
scribe the changes introduced by the presence of two deuterons 
in the PES of NH3 - H2, published recently ( Maret et al.|2009| l . 
We present in section |3] the collisional equations for the ND2H - 
H2 system. The cross sections and the corresponding reaction rate 
coefficients of the ND2H - H2 system are given in section |4] We 
present our conclusions in Section[5] 



We have recently put forward a potential energy surface (PES) of 
NH3 in interaction with H2. This rigid rotor PES is reported in full 
details in |Maret et al.|p009) . Here we need to transform this PES 
from the triple H isotopologue of ammonia to the doubly deuter- 
ated one, ND2H. For a given /nfer-molecular geometry, the inter- 
action between the molecules of ammonia and Ht depends on the 
various electronic charges (electrons and nuclei) as well as of the 
i>i?ra-molecular geometries of ammonia and molecular hydrogen. 
For the two isotopologues NH3 and ND2H the charge structure is 
identical, as long as we remain in the Bom-Oppenheimer approx- 
imation. However, the geometries are different in two aspects: (/) 
the intra-molecular distances change slightly and (ii) the position 
of the centre of mass of ND2H is shifted with respect to the centre 
of mass of NH3 and the principal axes of inertia are rotated. Let us 
treat both points in succession. 

The original PES for NH3-H2 was computed with both 
monomer geometries taken at their average value in the ground 
vibrational state ( [Maret et al.|2009T l. This has been demonstrated 
to be a much better approximation than taking equilibrium geome- 
tries: in the case of the similar H2O-H2 system, ' Faure et af1 ( |2005| l; 
[Valiron et al.| ( |200^ have found that employing state-averaged ge- 
ometries is a very good approximation for including vibrational ef- 
fects within a rigid rotor PES. We note that the importance of the 
zero-point vibrational correction was discussed earlier by |Meuw"ly| 
|& Bemish|p997) , in the N2H+ - He system. 

In particular the r(HH) distance for molecular hydrogen was 
shown to provide the largest vibrational correction. It was taken 
by [Maret et al.| ( |2009l l at its average value of 1.4488 Bohr. Sim- 
ilarly, the NH3 molecule was taken at its ground- state average 
geometry with the following parameter: r(NH) = 1.9512 Bohr, 
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Figure 1. Scheme of the coordinate transformation between the frame 
O'x'y'z' of ND2H, and the original frame Oxyz, of the NH3-H2 system, 
see Eqs. 1(2), j3). Only the Oxz and O'x'z' axes are represented, since the y 
axes do not change. Distances are not to scale. 



HNH = 107.38 deg. For ND2H, the equilibrium geometry is iden- 
tical to NH3 since the equilibrium bond lengths and angles depend 
only on the electronic structure, within the Bom-Oppenheimer ap- 
proximation. However the anharmonic averaging over the ground 
state vibrational state changes slightly those values, because of 
mass effects on the vibrational functions. For H2O/D2O, these ef- 
fects have been shown to be essentially negligible with respect to 
the shift of the center of mass jScribano et aL|20ro) . As a result, the 
above ground-state averaged geometry of NH3 can be employed for 
any isotopologue of ammonia and was adopted here. 

The second point deals with the fitting of the ab initio points 
onto a functional suited for the scattering calculations. The PES 
points resulting from the ab initio procedure are described in the 
following coordinate set: R, the distance between the centre of mass 
of NH3 and the centre of mass of H2, 9, if>, the corresponding spher- 
ical angles. For NH3, d is defined with respect to Oz, the C3,, axis 
of symmetry of NH3, with the nitrogen atom at positive values of z- 
The angle (f> is the polar coordinate in the perpendicular Oxy plane, 
with one of the hydrogen atoms of NH3 lying in the Ozx plane. Sim- 
ilarly, the H2 molecule is oriented by the spherical angles du, 0h, 
in a reference frame OnXnynZn, parallel to Oxyz , with its origin 
at the centre of mass of the H2 molecule (see Fig. 1). In order to 
define the new coordinate system for NDiH, we must first anchor 
the new centre of mass in the old coordinate system, then perform 
the relevant rotation of the Oxz axes around the invariant Oy axis, 
which remains perpendicular to the plane of symmetry of the ND2H 
molecule. If O'x'y'z' is the new spherical frame centered on the O' 
centre of mass on ND2H, and 7 the rotation of the O'z' axis around 
the O'y' axis (O'y' is parallel to Oy), we have according to [Cohen] 
|& Pickett| ( [T9g2l i : 



xo' = -0.095957 Bohr 
zo' = -0.062152 Bohr 
7 = 11.30 degrees 



(1) 



It is enough to translate the centre of mass and perform the 
relevant rotation 7 around Oy and OhVh to get the new coordinate 
system in which the ND2H-H2 PES should be expressed (see Fig.[TJ 
. The change of coordinates is given by: 



X = x' cos 7 - z' sin 7 -I- Xq' 

y=y' 

z = x' sin 7 + z' cos 7 + zo' 



(2) 



that O'x'y'z' and OnxJjyJjZjj remain parallel and Ojj is now located 
with respect to O'x'y'z' at a distance R'. This corresponds to a pure 
rotation of 7 along the yu axis: 



X[j = Xfj cos y — Zjj sin 7 

>'«=>'« 
Zh = x'fj sin y + z'u cos 7 



(3) 



Eventually, the functional form we take is identical to the one 
for water-hydrogen ^Phillips et al.|1994[[Wiron et al.|2008) . The 
spherical angles, all primed, refer to the new frames O'x'y'z' and 

V {R',e',</)',e'^^, cp'ii) = ^ v,,,,,,„(i?')?;,/,/,„ <P', <P'h) • (4) 



The H2 internal coordinates are transformed in a similar way, so 



where the function are explicitly given in Phillips et al.| 

(19941; WTiron et al. ("50081. As previously ( [Valiron et al.||2008[ 
|Maret et al.„2009j , we selected iteratively all statistically signifi- 
cant terms vi^i^i^(R') using the procedure described in | Valiron et al.| 
( |2008^ . The final expansion includes anisotropics up to ;i = 10 for 
ND2H and /2=4 for H2, resulting in 197 J/j/j/,,, functions. 



3 COLLISIONAL TREATMENT 

In the low temperature environments where ND2H has been de- 
tected the kinetic temperature is below 20 K ( [Roueff et al.|2000l > 
and the main constituent, H2, is expected to lie in its ground ro- 
tational (para) state, 7 = ( [Troscompt et al.|2009) . The calcula- 
tions presented below were restricted to collisions between NDiH 
and para-H2(/ = 0), i.e. by neglecting the J = 2 (closed) channel. 
This latter was actually found to change the cross-sections, at very 
low energies, by up to a factor of 3, as observed for NH3-H2 (see 
references in [Maret et al.[ ( (2009| l). Its effect on the average cross- 
sections and rate coefficients is however much smaller, typically 
30 %, which is therefore the typical accuracy of the rates below. 
The rotational constants of ND2H were taken from [Coudert et aE] 
l [1986| . The resulting rotational energy levels in are given Table [T] 
The reduced coUisional mass is /j = 1.822684 amu. 

Full quantum close coupling (CC) scattering calculations were 
performed with help of the Molscat program ( [Hutson & Green] 
[1995[ l which computes scattering matrices (5 matrices) and com- 
bines these in order to determine elastic and inelastic cross sec- 
tions, at defined total energies E. The convergence of the inelas- 
tic cross sections was checked to be better than 5%, for the dif- 
ferent rotational transitions investigated here. The different col- 
lision parameters defined in Molscat are given the default ones. 
We used the following values : RMIN = 3.0, INTFLG = 6 , 
STEPS = 10.0 - 30.0. Calculations have been performed for to- 
tal energies between 10 cm"' and 430 cm"'. The energy step has 
been varied with increasing collision energy. It was 0.25 cm"' for 
a total energy between 10 cm"' and 40 cm"', 1 cm"' between 40 
cm"' and 80 cm"', 5 cm"' between 80 cm"' and 100 cm"' and 20 
cm"' above 100 cm"' . The size of the rotational basis set of ND2H 
was also varied with energy. For E < 100 cm"', y(ND2H) < 10; 
at £ = 120cm"', y(ND2H) < 12. For the CS calculations at 
E > 150 cm"', i(ND2H) < 16. 

In order to limit the time of computation needed for a scatter- 
ing calculation, some customary approximations of the collisional 
treatment are implemented in the Molscat program such as the Cou- 
pled State (CS) approximation introduced by | McGuire & KourT] 
( [1974^ . In this approximation, scattering equations are written in 
the body-fixed frame which is rotating. Then the projection of the 
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Table 1. Rotational level energies for ND2H, as used in our computation. 





Energy (cm 


Ooo 





lot 


9.09393 


In 


11.19743 


lio 


12.78456 


2o2 


26.65954 


2l2 


27.79543 


2ii 


32.55350 


221 


38.86400 


220 


39.48154 



orbital angular momentum / is restricted to the value m; = 0. The 
accuracy of this method is questionable but it was shown to agree 
with CC calculations within typically ~ 30% for the present sys- 
tem, as shown below. Here we resorted to the CS approximation 
for total energies E > 150 cm"'. As a result, at the low tempera- 
tures at which the rates were calculated (T < 35 K), the sensitivity 
to the CS calculations was checked to be marginal. 



4 RESULTS 

Figure[2]displays the collisional de-excitation cross sections corre- 
sponding to the transition in — » Ooo and lio — » loi as a function 
of the relative kinetic energy between ND2H and H2, both at the 
CC and CS level. These transitions are observable from the ground 
at 335 and 110 GHz and are potentially important probes of cold 
prestellar cores ( |Lis et al.|[2006[ |Gerin et ar]|2006[ [Roueff et aT] 
|2000 l). For comparison, also cross sections with He computed at 
the CS level, from Machin & Roueff ( 2007j , are shown. A com- 
plex resonance structure is seen for energies smaller than about 100 
cm"' with broad and narrow features, resulting in large variations 
of the cross sections. The sharp maxima are due to the opening of 
new collision channels and correspond to both Feshbach and shape 
resonances. This is in contrast with the He-NHD2 collision case, 
with a much shallower potential well. Even if the symmetry of the 
two problems are identical, differences arise in their cross-section 
structures and cross section values. A similar effect was found re- 
cently at much higher energies, for the differential inelastic cross 
section of water-H2, 7 = compared with water-He |Yang|2010j . 
As is always been observed, the resonances quickly disappear with 
increasing energy. The cross sections relative to Helium are sig- 
nificantly smaller in the energy range displayed in Figure |2] How- 
ever, they converge to those involving para-H2 at collision energies 
above ~80 cm"'. 

Knowing the importance of the hyperfine quadrupolar effect 
for '''N, nuclear spin / = 1, extending those cross section calcula- 
tions to include explicitely the hyperfine structure could be of inter- 
est. Hyperfine structure has indeed been observed with some preci- 
sion in the lower transitions I jGerin et al.|2006[ l. However, because 
of the low symmetry of NHD2 (asymmetric top), no close formula 
exist in the literature that could connect state to state hyperfine rates 
to state to state global rotational rates. One should resort to a full 
sumation of S matrix elements in order to get all detailed hyperfine 
rates. Also, because of the small column density of NHD2, the dif- 
ferential effects on the various hyperfine lines should remain small 
Paniel et al.|2006t . 



The collisional (de)excitation rates Rasa, function of temper- 
ature T were obtained by a Maxwell-Boltzmann averaging of the 
cross sections times the relative velocity. In the equation below, 
Eti„ . is the relative kinetic energy: 



SkeT 

Tip. 



(5) 



kg is the Boltzmann constant and ji is the reduced mass of the ND2H 
- H2 system. Cross sections having been computed for collision en- 
ergies up to above 390 cm"', the rate coefficients may be calcu- 
lated very safely for temperatures between 5 K and 35 K which are 
relevant for cold prestellar cores. Table |2] gives the de-excitation 
rate coefficients for the 9 first rotational levels as well as the fitting 
parameters allowing to display the rate coefficients for any tem- 
perature between 5 and 35 K via the following expression l |Faure,| 
Gorfinkiel, & Tennyson|2004|l, with R in 10" ' 



n/6 



(6) 



All fits have been verified to be precise to better than 1 %. It should 
be noted that the following inelastic transitions are strictly forbid- 
den in the CS approximation (owing to spurious selection rules): 
loi — > Ooo J 2i2 — » Ooo and 2ii — > Ooo- We see that the corre- 
sponding CC rates are not negligible, although about one order of 
magnitude smaller than the largest rates. There is therefore no ob- 
vious collisional propensity rules for the presented transitions. We 
note that since the dipole moment has two non-zero projections (//(, 
and /Jc), the corresponding radiative selection rules are AJ= 0, 1 and 
AK„ = ± 1, AKc = ±1 for transitions along fii, and AKa =±1, AK^ = 
0, ±2 for transitions along fic- The Einstein coefficients have been 
computed by Coudert and Roueff] ( |2006| ) and recalled in |Machin &| 



|Roueff| ( |2007 



5 DISCUSSION AND CONCLUSIONS 

The comparison of the present rates with the values obtained by 
[Machin & Roueff] ( |2007^ where Helium is the collider is informa- 
tive. The collision rates obtained in the present work when para-H2 
(J=0) is the perturber are typically one order of magnitude larger 
than those obtained with Helium. The origin of this discrepancy 
lies in the significantly different potential wells of the NH3-H2 and 
NHj-He PES, as discussed by jMaret et ar]p009) . The usual scal- 
ing between the collision rates obtained from the ratio of the re- 
duced masses is clearly not accurate and the calculations have to be 
performed explicitly in order to obtain the appropriate rates with 
molecular hydrogen. We display in Table [3] the resulting critical 
densities for the para and ortho transitions. These were obtained 
by simply dividing the Einstein coefficient of a given transition by 
the corresponding de-excitation rate coefficient at 10 K. Note that 
other definitions are possible ( [Maret et aH2 009 1. These values set 
the limit of the perturber density above which the transitions be- 
come thermalized. We see that the critical densities corresponding 
to the two transitions In — > Ooo and lio — > Ooo are different by 
about one order of magnitude at 10 K, which is about a factor of 
3 less than the ratio of the Einstein coefficients. Such differences 
may have profound consequences on the interpretation of these two 
transitions which are both observable from the ground ( jLis et alT] 
[2006] [Gerin et al.|2006) at 335 and 389 GHz respectively, as the 
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Figure 2. Collisional de-excitation cross sections as a function of tlie relative kinetic energy. Left: transition In — » Ooo of ND2H witli Helium in CS 
approximation (red) ((Machi n & Roueff|2007) ) and para H2 (present work) in CC (blue) and CS (green) approximations. Right : transition lio — » loi of 
ND2H with Helium in CS approximation (red)( jMachin & Roueff|2007) l) and para H2 (present work) in CC (blue) and CS (green) approximations 



Table 3. Calculated critical densities for the first ortho and para transitions 
ofND2H withH2 at 10 K. 



Transition Critical density (cm ^) 







ortho 


para 


In 


-^Ooo 


1.17 10* 


1.34 10* 


111 


^ loi 


3.48 10-* 


6.63 lO"* 


lio 


^Ooo 


1.10 10^ 


1.01 10' 


lio 


^ loi 


5.19 m* 


5.94 lO* 


2o2 


^ In 


6.34 lO'' 


6.76 10^ 


2o2 


^ lio 


5.33 10'' 


5.66 10* 


2l2 


-> 2o2 


2.31 lO* 


9.26 10^^ 


2ii 


^ loi 


1.14 10** 


1.18 10** 


2ii 


^ lio 


1.33 10^ 


1.45 10' 


2ii 


— > 2o2 


7.59 lO"* 


7.12 lO* 


221 


-» 111 


1.80 10** 


1.79 10** 


221 


-> lio 


4.17 10'' 


1.77 10' 


221 


^2i2 


2.70 10"^ 


1.94 10^ 


221 


^2ii 


1.39 10'' 


1.23 10* 


220 


-> 111 


5.78 10^ 


2.26 10' 


220 


-» lio 


1.47 10** 


1.49 10** 


220 


^2i2 


5.65 10'' 


6.44 10* 


220 


-*2ii 


3.96 lO'' 


2.48 10^ 



corresponding upper levels may not be accounted for by a single 
excitation temperature. 
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Table 2. Deexcitation rate coefficients for ND2H-para H2 (J = 0) collisions. Fits to formula j6| are given in the last 4 columns. 









Rate coefficients R,y 


( 10-" cm' 








Fit coefficients 




Temperature 


5K 


lOK 


15K 


20K 


25K 


30K 


35K 


a„ 


1 


2 


3 


Iqi 


Ooo 


0.572 


0.436 


0.351 


0.29 


0.245 


0.211 


0.183 


19.71 


-46.57 


37.81 


-10.74 


In 


—> Ooo 


1 14 


1.10 


1.03 


0.938 


854 


78 1 


717 




-51.13 


38.73 


-9 74 


111 


^ loi 


7 S7 
/ .3 / 


6.79 


6.14 


5.62 


S 7 1 






u.uu 


-3.51 


5.63 


1 

- i . JO 


lio 


Ooo 




4 39 


4 05 






^ 90 


J. Lj 


1 89 


- / . JO 


8 07 


-2 04 


lio 


— > loi 


1.49 


1.34 


1.23 


1.12 


1.03 


0.953 


0.887 


15.74 


-35.56 


27.24 


-6.91 


lio 


^ 111 


1.37 


1.25 


1.14 


1.05 


0.959 


0.887 


0.826 


14.76 


-33.90 


26.31 


-6.76 


2o2 


Ooo 


824 


0.782 


0.71 


0.64 


581 


532 


493 


O.Jo 


-19.33 


17.89 




2o2 


^ loi 




0.842 


0.709 


0.613 


537 


476 


425 




-63.93 


46.67 


-11 72 


2o2 


^ In 


J. 1 7 


3.06 


2.97 


2.90 


Z.OJ 


9 77 
z. / / 


9 11 


7 '>7 


-5.88 


4.76 


C\ 89 
-U.oZ 




^ lio 


1 9Q 


1.12 


1.00 


0.904 


821 


752 


694 


1 

LJ.yo 


-36.13 


27.88 


7 9'^ 

/ .Z J 


2l2 


^ Ooo 


0.275 


0.272 


0.234 


0.199 


0.171 


0.148 


0.13 


6.25 


-24.01 


24.97 


-8.41 


2l2 


-> loi 


j.y J 


3.84 


3.73 


3.61 


3 52 


3 44 


J.J 1 


1 09 


-3.63 


3.67 


-U.Jo 


2l2 


^ In 


808 


0.868 


0.823 


0.753 


684 


623 


569 


15 88 


-40.31 


32.97 


-8 84 


2l2 


lio 


n 7^Q 


0.787 


0.744 


0.683 


621 


565 


W.J 10 


17 41 


-43.09 


34.66 


Q 99 


2l2 


^ 2o2 




3.77 


3.97 


4.03 


4 03 


4 01 


'\ QQ 
J.yy 




-25.60 


17.74 


-3 40 


2ii 


^ Ooo 


104 


0.116 


0.102 


0.0875 


0749 


0647 


0564 


2 64 


-19.41 


23.36 


-8 68 


2ii 


-> loi 


1.95 


2.22 


2.36 


2.40 


2.40 


2.38 


2.35 


16.14 


-34.63 


24.05 


-5.07 


2ii 


^ 111 


3.72 


4.02 


4.14 


4.14 


4.08 


3.98 


3.88 


18.59 


-39.65 


27.74 


-5.77 


2ii 


lio 


464 


0.489 


0.464 


0.428 


\J.jyj 


362 


\).Jj\) 


10 75 


-28.86 


24.57 


-7 06 


2ii 




2 78 


2.99 


3.00 


2.96 


2 90 


2 85 


2 81 


1 09 
i .wz 


-0.82 


2.40 


-0 45 


2ii 


^ 2i2 


1 08 


1.25 


1.23 


1.15 


1 07 


1 00 


\).y^J 


Z.JO 


-13.06 


14.38 


-4 38 


221 


^Ooo 
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